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q \ ABSTRACT 

Aims. We use data from observational cosmology to put constraints on higher- 

. „. 

dimensional extensions of general relativity in which the effective four-dimensional 

dark-energy density (or cosmological "constant") decays with time. 

Methods. In particular we study the implications of this decaying dark energy for the 

OO ' 

age of the universe, large-scale structure formation, big-bang nucleosynthesis and the 



^1- 

cn 
o 



X 



magnitude-redshift relation for Type la supernovae. 

Results. Two of these tests (age and the magnitude-redshift relation) place modest lower 
limits on the free parameter of the theory, a cosmological length scale L akin to the 



de Sitter radius. These limits will improve if experimental uncertainties on supernova 
magnitudes can be reduced around z ~ 1. 
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1. Introduction 

In standard general relativity, dark energy is interchangeable with Einstein's cosmological 
constant A and the dark-energy density p A = A/(8irG) is constant. Observation tells us 
that, at the present time, this latter quantity is nonzero but many orders of magnitude 
smaller than expected based on calculations in quantum field theory. This mismatch has 
led some theorists to look at alternatives to standard general relativity in which A is 
a dynamical parameter whose value is not constant but might have decayed from large 
values in the early universe to those we see today (Overduin & Cooperstock 1998). 

Here we look at one such alternative, the minimal extension of Einstein's theory to 
manifolds with one additional noncompact spacelike dimension (Overduin & Wesson 1997). 
Higher-dimensional gravity has been shown to be compatible with solar-system tests in- 
cluding the classical tests of general relativity (Kalligas et al. 1995, Liu & Overduin 2000) 
and experimental limits on violations of the equivalence principle (Overduin 2000). We wish 
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to check here whether it is also consistent with basic cosmological tests that are sensitive 
to the A term. An exhaustive treatment of all cosmological data is beyond the scope of 
this introductory study. 

Mashhoon & Wesson (2004) have shown that a gauge transformation in higher- 
dimensional gravity converts A from a constant of nature to a gauge-dependent measure 
of four-dimensional vacuum-energy density. Under reasonable assumptions (conformal flat- 
ness and geodesic motion in five dimensions), A decays exponentially with cosmic time. We 
investigate the cosmological implications of this kind of vacuum decay by solving numer- 
ically for the age t(z) of the universe as a function of redshift. The bulk of the vacuum 
decays takes place near z ~ 1, too late to affect big-bang nucleosynthesis (BBN) or large- 
scale structure formation (LSS). Lower limits on the age of the universe, however, put weak 
constraints on the primary free parameter of the theory (a de Sitter-like length parameter 
L), and these bounds are improved somewhat by data on the magnitude-redshift relation 
for Type la supernovae. Taken together with the solar-system tests, we tentatively conclude 
that available astrophysical data are consistent with a universe with one (or more) extra 
dimensions. 



2. Higher-Dimensional Cosmology 

A starting point for cosmological investigations in 5D theory is the metric in canonical 
form (Mashhoon et al. 1994): 

£ 2 

dS 2 = — [g^{x x ,t)dx^dx u ] - dl 2 , (1) 
±j 

where L is a constant with dimensions of length (akin to the de Sitter radius Ld — -\/3/A 
in standard cosmology). The 5D line element contains the 4D one: 

ds 2 = g flv (x x ,£)dx fJ 'dx v . (2) 

There is no loss of generality to this point; hve available degrees of coordinate freedom 
have been used to set the electromagnetic potentials (#4^) to zero and to set the scalar 
potential ((744) to a constant in Eq. ([T]). It is, however, necessary to retain ^-dependence in 
the 4D metric tensor in order to preserve this generality (Overduin & Wesson 1997). 

Under the restriction to 5D conformal flatness, and the natural assumption that all 
test particles (massive as well as massless) move along null geodesies in 5D (i.e., dS 2 = 0), 
Mashhoon & Wesson (2004) have shown that A in 4D drops exponentially with proper 
time s: 

3 1 

1? (l-e=F«/£)2 ' ( ) 

Physically, this variation arises because we require the 5D field equations to satisfy general 
covariance in five, not four dimensions. The canonical metric ([1} is invariant with respect 
to translations along the ^-axis, so t- or gauge-dependence then necessarily appears in the 
4D field equations. We have used the 4D metric to re-express this dependence in terms of 
proper time s rather than I. There are two cases: in the first (' — ' sign in the exponent), 
A decays asymptotically to the small finite value 3/L 2 as s — > 00, while in the second ('+' 
sign) it vanishes in this limit. Measurements tell us that A is small at present, but are not 
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precise enough to discriminate between a constant value and one that is still decaying on 
cosmological timescales. Therefore we retain both possibilities in what follows. 

One way to constrain proposals of this kind is to ask what A decays into. If matter or 
radiation, then strong constraints can be placed on the theory using experimental limits 
on cosmic background radiation (Overduin et al. 1993, Overduin & Wesson 2003, 2004). In 
this paper, we will investigate the more conservative scenario in which A docs not decay 
into matter or radiation, so that the evolution of the matter density (p M ) and radiation 
energy density (p R ) proceed just as in standard 4D cosmology. It should be noted that such 
an assumption requires, in principle, the existence of some other field to which A is coupled 
and into which its energy can be transferred, in accordance with the 4D conservation law 
(Overduin & Cooperstock 1998): 

V"(8.GV-A V )=0. (4) 

In order to compare Eq. ([3]) to observation, we need to convert from proper time s to 
ordinary cosmic time t. A physically-motivated argument due to Mashhoon is the following: 
translations along the £-axis do not only introduce gauge-dependence into A; they also give 
rise to an apparent "fifth force" in the equations of motion (Mashhoon & Wesson 2004, 
Wesson 2005). This force can however be made to vanish by an appropriate choice of affinc 
parameter, so making the motion geodesic in the usual 4D sense (Seahra & Wesson 2001). 
This choice can be shown to lead to the following relation between s and t, assuming 
that the motion in 5D is null, that the galaxies are comoving in 4D spacetime, and that 
dt/ds > over < s < oo: 

e ±s / L = l±{l/a)e t / L . (5) 

Here a (expressed this way for later convenience) is an unknown positive constant whose 
value is to be fixed by cosmological boundary conditions, and the signs are to be read in 
the same way as in Eq. Q; i.e., the upper ('+') sign here corresponds to the upper (' — ') 
sign there. Putting Eq. © into the decay law fl3|), we find that: 

A(t) = -^(liae-^) 2 , (6) 

where the sign order again corresponds to that in Eqs. and ([SJ> . 

Eq. © provides the starting-point for our investigation. Motivated by the 5D results, 
we study dark-energy decay of the functional form given by Eq. ^ in the context of 4D 
cosmology from t = at the big bang to t — > oo. This kind of exponential dark-energy decay 
appears to be unique in the literature (Overduin & Cooperstock 1998). In the limit t — > 0, 
we observe that A originates with a finite value of 3(1 ± a) 2 /L 2 and decays asymptotically 
toward 3/L 2 as t — > oo. Motivated by the cosmological-constant problem, we might expect 
very large values of a with steep drop-offs at early times. We shall find, however, that a 
is typically within a few orders of magnitude of unity, rendering departures from standard 
4D cosmology rather mild. 
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3. Decaying Dark Energy 

Evaluating Eq. © at the present time t , restoring physical units and expressing the results 
in terms of the critical density, we obtain: 

where H is the present value of Hubble's parameter. It will be convenient to rescale all 
dynamical parameters in dimensionless terms, so wc define: 

T = H a t , C = H L/c , H = H/H . (8) 

In terms of these quantities Eq. ([7]) determines the age of the universe r = H 1 in terms 
of the two free parameters a and C as follows: 

We note from Eq. ([9]) that the two cases corresponding to the '+' and ' — ' signs are sep- 



arated by a "critical" case with £ crit = 1/ \/^a,o- F° r values of C > C crit we must use 
the '+' solution, while the ' — ' solution is operative if L < C am . [An alternative solution, 
corresponding to the negative root of Eq. ([7]), is also available in principle, since Eq. © 
has three possible square roots. However, the root that results in the alternative solution 
should be excluded as it does not have the same limit as Eq. (J6j) for t — » oo. This explains 
why the alternative solution is not taken into consideration here.] In the limit C — > C otH 
precisely, it is apparent from Eqs. ([7|) with the definitions ([8|) that a = and the theory 
goes over to standard cosmology with A =const. The theory also goes over to standard 
cosmology in the limit L — > oo, which from Eq. ([7]) means that a must go as L 2 for large 
L, as we will confirm numerically below. The physical length corresponding to £ crit is just 



the de Sitter radius of standard cosmology, L cm = c/{H ^/Q. A ), which takes the value 
L crit = 4.9 Gpc for WMAP values of H and f) AiD (Spergel et al. 2003). 

We now fix the value of a by requiring consistency between Eq. (J9j) and the age of the 
universe r as obtained by numerical integration of the Friedmann-Lemaitre equation: 

T{z)=n ~L (l + Z >)H [z >r^)}- (10) 
Here z is redshift and the standard expression for Hubble's parameter is modified following 
Eqs. © and © so that: 



H[z,t(z)} 



n Ml0 (i + z? + (i - ^M, - Ov,o)(l + zf 



C- 2 l±ae- T ^/ c 



1/2 



(11) 



Consistency is to be enforced by putting the boundary condition r(oo) = into Eq. (flT))) : 
i.e., for each value of L we solve numerically for a value of a satisfying: 

dz ' c\J \ . (12) 



/ (i + z')h[z',t(z')] \c^nz-i t 

We adopt WMAP values for the present values of the density parameters, assuming spatial 
flatness: Q M , = 0.135/^ = 0.27 (with h = HjlOQ km s" 1 Mpc" 1 = 0.71) and fl A , = 



J.M. Overduin et al.: Decaying Dark Energy in Higher-Dimensional Gravity 5 

1 — fi MiD = 0.73 (Spergel et al. 2003). That leaves us with only one adjustable parameter 
in the theory: the de Sitter- like length L. 

Fig. 1 (top) shows how this cosmological consistency requirement produces values of a 
between approximately 0.1 and 1000 for a range of L- values with L > 5 Gpc. As expected, 
a — > in the limit L — > L crit and a oc L 2 for large L. 

Fig. 1 (bottom) shows the corresponding age of the universe t = t /H [as obtained 
from either Eq. © or Eq. ([TUfl in the limit z — ► oo] for the same range of values of L as 
in Fig. 1 (top). The shape of this plot can be understood physically as follows: The age 
of the universe goes over to that of a standard flat ACDM model with fi M ,o — 0-27 (i.e., 
13.7 Gyr) in both the "critical" case L — L CIit and the limit L —> oo, as expected. To 
the right and left of the critical case are the '+' and ' — ' regimes in Eq. (jSJ) respectively. 
The presence of a local minimum and a maximum here is due to the L-dependence of the 
dark-energy density: the behavior near these extrema is dominated by the exponential term 
±aexp(— ct/L) in Eq. ([6]), whose magnitude in this region is comparable to unity. 

To understand why the age of the universe in this theory is always lower than that 
of standard cosmology for L > L crit , whereas it can climb somewhat above this value for 
L < L„ it , we recall that the dark-energy density is always pinned at its observed value at 
present. In the '+' regime (L > L crit ), it can only grow in the past direction relative to 
its value in standard cosmology, increasing the total density of matter plus dark energy. 
Since the square of the expansion rate is proportional to total density (from the Friedmann 
equation) the expansion rate goes up and the expansion timescale goes down. Hence, for 
given values of £l M ,a and ^a,o the age of the universe in this theory is always smaller than 
that in standard cosmology for L > L crit . In the ' — ' regime, by contrast, there is a range of 
L- values for which the dark-energy density in the past direction is slightly lower than that 
in the equivalent standard cosmology (i.e., with the same boundary conditions); hence one 
obtains a slightly older universe in this region. 

These departures from standard cosmology can be significant: for L > L crit the age of 
the universe drops from 13.7 to a local minimum of 13.0 Gyr at L mia — 9.2 Gpc, while for 
L < L crit it climbs to a possible maximum of 16.9 Gyr at L max = 3.5 Gpc, before dropping 
rapidly toward zero as L — > (assuming WMAP values of f2 M . and O Ai0 as usual). This 
curve is sufficient to set a robust lower bound L > 2.2 Gpc from the fact that the age of the 
oldest observed stars in globular clusters are at least 11 Gyr (Wanajo et al. 2002, Schatz 
et al. 2002). The range of acceptable L values could be further constrained by other limits 
on the age of the universe from a variety of observations, including upper bounds (though 
the latter are necessarily less robust since we may not be able to see the oldest members 
of any target population). In any case we will obtain stronger constraints below. 

As a first step we use Eq. (|10p to compute the age t(z) = t(z)/H as a function of 
redshift. The presence of explicit time-dependence under the integral sign makes this a 
difficult equation analytically, but it can be solved numerically. Results are shown in Fig. 2 
(top). The L = L cm ,L — L min and L = L max cases are again indicated together with a 
representative sampling of other models. 

The quantity t(z) plotted in Fig. 2 (top) is helpful in checking whether the theory 
is consistent with the formation of LSS by gravitational instability in the early matter- 
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dominated era and BBN at the end of the radiation-dominated era. For these processes to 
work successfully, the density of dark energy must not rise so quickly in the past direction 
that it becomes comparable to or greater than the matter density during LSS formation, 
or comparable to or greater than the radiation energy density during BBN. These latter 
two quantities are given by: 

Pn(z) = p„it,o ^r,o(1 + z) 4 



where p crit ,o = 3Hq/(8ttG) = 1.88 x 10~ 29 h'g g cm -3 is the critical density and fi Ri0 = 
4.17 x 10 _5 //ig from COBE data plus standard neutrino physics (Overduin & Wesson 
2003). By comparison, the density of dark energy is given by Eq. ^ as: 



with t(z) as shown in Fig. 2 (top). 

Inserting t(z) into the dark-energy density (|14[) , and comparing with the matter and 
radiation energy densities, we obtain the plots shown in Fig. 2 (bottom). A representative 
sampling of models with L > 3.5 Gpc is included, with the L = L crit , L mln and L max cases 
labelled as before. It is clear from these plots that departures from standard cosmology are 
too small to have a significant effect on either LSS or BBN. For L < L crit the density of 
dark energy drops relative to that in standard cosmology. For L > L crit , even the largest 
possible increase in dark-energy density (corresponding to the shortest-lived universe with 
L = L min ) is far too modest relative to the matter density at z > 1 to interfere with structure 
formation, and completely negligible relative to radiation-energy density at z > 1000. Thus 
these tests do not place meaningful constraints on the theory. 

Fig. 2 (bottom) shows that the largest departures from standard theory are found near 
z ~ 1, raising the possibility that stronger constraints might be obtained by use of the 
SNIa magnitude-redshift relation. Supernovae are now being routinely monitored at z ~ 1 
(Riess et al. 2006), providing a sensitive testbed for alternative theories of gravity with 
time- varying dark-energy density (see for example Fukui 2006). 

The magnitude-redshift formulae are derived in the Appendix. We focus on the mag- 
nitude residual Am(z), or difference in apparent magnitude relative to a fiducial model, 
which we take here as the standard flat ACDM model with WMAP values of fi Mi0 and 
r2 A>0 . Predictions are plotted as curves (for various values of L) in Fig. 3, where they are 
compared with measurements for 92 medium-redshift SNIa at z > 0.1 by Tonry et al. 
(2003) and 23 high-redshift SNIa at z ~ 1 as compiled by Riess et al. (2006). 

The heavy solid line in Fig. 3 indicates the fiducial or standard ACDM model (straight 
line Am = 0), which overlaps with the present theory in the case L = L CIit and all cases 
with L > 300 Gpc (to within the precision of the plot). For smaller values of L, the theory 
begins to depart from standard cosmology with maximum deviations near z ~ 1, confirming 
the usefulness of SNIa as probes of the theory. For L near (or greater than) the critical 
value i crit = 4.9 Gpc, the experimental uncertainties are too large to discriminate usefully 
between theoretical values of L. The data are, however, good enough to disfavor smaller 
values of L, improving significantly on the age constraint and tightening observational 
bounds on the theory to L > 4.3 Gpc. 



Pm(z) = Pcrit,O^M,o(l + Z) 3 , 



(13) 




(14) 
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4. Discussion and Conclusions 

We have taken the basic extension of general relativity from 4 to 5 dimensions and asked 
what observational consequences follow from its decaying cosmological "constant," or den- 
sity of dark energy. The four main consequences involve the age of the universe, structure 
formation, nucleosynthesis and the magnitude-redshift relation for Type la supernovae. 
There are, of course, large literatures on all four of these subjects. We have therefore pre- 
sented our results as possible departures from the current standard model. The theory 
is consistent with all four classes of data at present. The best way to separate 4D and 
5D cosmology by observational means in the future would appear to be by use of better 
supernovae data at z > 1. 

Another possible test of the theory might come from analysis of the angular power 
spectrum of fluctuations in the cosmic microwave background (CMB). Qualitative con- 
siderations, however, suggest that the sensitivity of such a test would not be competitive 
with those discussed above for the kind of theory considered here. The main feature of 
the CMB power spectrum is the angular position of the first acoustic peak, £ poak . This 
quantity depends only weakly on dark-energy density [see Fig. 1 of White (1998) or the 
analytic approximation in Cornish (2000)] . What £ poak really measures is the sum of mat- 
ter and dark-energy densities, i.e. spatial curvature. We have assumed throughout that 
O Mj0 + A|C = 1 (i.e. k = 0), as in the standard ACDM model, so we would not expect 
a significant shift in £ poak . Fig. 2 (bottom) shows that dark-energy density changes by at 
most ~ ±40% between zwl and z ~ 10 relative to standard ACDM cosmology. Following 
Cornish (2000) the change in A might shift £ poak by at most ~ 40%/35 ~ 1%, compara- 
ble to the current level of experimental uncertainty in this parameter (Page et al. 2003). 
Physically, the reason why the CMB constrains this theory less strongly than SNIa is be- 
cause it probes higher rcdshifts where the density of matter is so much higher than that 
of dark energy that a modest change in the latter has little effect. More detailed study is 
warranted, however, particularly in light of the increase in experimental precision expected 
from the Planck mission. We hope to return to this issue in future work. 

These cosmological results are complementary to earlier ones based on the classical 
solar-system tests and ones involving the equivalence principle. It has been known for 
some time that the basic 5D extension of 4D Einstein theory is consistent with these tests 
(Kalligas et al. 1995, Liu & Overduin 2000, Overduin 2000; Wesson 2006 for review). It is 
important to recall here that standard cosmological models which are curved in 4D may 
be smoothly embedded in models which are flat in 5D. [For a review see the books by 
Wesson (1999,2006); for an account of the embeddings from an astrophysical viewpoint see 
Lachieze-Rey (2000).] The data may therefore be suggesting not only that the universe has 
a fifth dimension, but that its structure may be much simpler than previously thought. 

Appendix 

The apparent magnitude m of a source at redshift z is defined by: 



m(z) = M + K(z) + 51og[d L (z)/10pc\ , 



(15) 
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where M is absolute magnitude, K{z) is the K-correction due to frequency shift and the 
luminosity distance d L (z) is given by: 

r dz' ' 
Jo n(z')_ • 

Here the constant ak and function Sk are defined so that ak 



c{l + z) 

d^z) = — S k 

n <?k 



(16) 



{\J$Im,o — ^a,o — 1 5 lj \/l — ^m,o — ^a.o} and <Sfc[X] = {sinX, X, sinhX} respectively 
for k = {+1, 0, —1}. Eq. (fl~5|) contains terms such as M, K(z) and H Q that are independent 
of the background cosmology and hence not of interest to us here. We remove those terms 
by focusing, not on the apparent magnitude itself, but on the difference or "residual" 
magnitude relative to a fiducial model, which we take here to be the standard flat ACDM 
model with WMAP values of f2 M0 and f2 A , . That is, we focus on the observational 
quantity 



Am(z) = m(z) — m m (z) = 5 log 



V{z) 



V tid {z) 



(17) 



where 
V{z) 



dz 1 



H[z',t(z')} 



r - dzi as) 

Jo v/O m , (1+^) 3 + ^a,o 
and where the modified Hubble parameter H is specified as before by Eq. (fTT|) . To compare 
our predictions with the magnitude residuals measured by Tonry et al. (2003) for 92 SNIa 
at z > 0.1, we write Am(z) = 5(y ± Sy) — 5 log[c(l + z)T> {id (z)], where z, y and Sy are read 
from columns 7, 8 and 9 of Table 15 in that paper [y = log(d L H )]. To incorporate the 
new and invaluable survey of 23 SNIa at z ~ 1 compiled by Rless et al. (2006), we note 
that y = log(d L H a ) = ^(/^o — C) where /i is read from column 3 of that paper and C is a 
constant whose value is fixed by requiring that both samples give consistent values of d L 
for SN1997ff at z = 1.755, implying that C = 15.825. 
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Fig. 1. (Top) Values of the constant a as a function of L, assuming WMAP values of f2 Mi0 
and Q A . (Bottom) Age of the universe as a function of L for the same values of fl M ,o and 
r2 A>0 . Labelled points correspond to the critical case separating the two classes of solutions 
of Eq. ([9]), and to universes of maximum and minimum age in the two regimes (see text 
for discussion) 
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Fig. 2. (Top) Age of the universe as a function of redshift for various values of L, assuming 
WMAP values of £1 M and Ai0 . (Bottom) Densities of radiation, matter and dark energy 
as a function of redshift for various values of L, assuming the COBE value of R|C and 
WMAP values of f2 Mj0 and O Aj0 . The three labelled points correspond to universes with 
critical, maximum and minimum values of L, as in Fig. 1 (bottom). 
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Fig. 3. Magnitude-redshift relation for various values of L, as compared to observational 
data on Type la supernovae 



